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Abstract. Extensive optical and infrared photometry as well as low and high resolution spectroscopy are used 
as inputs in deriving robust estimates of the reddening, distance and nature of the progenitor of V838 Mon, the 
2002 outbursting event that produced a most spectacular light-echo. The reddening affecting V838 Mon is found 
to obey the i?y=3.1 law and amounts to (i) Eb-v=0.86 from the interstellar Nal and KI lines, (ii) Eb-v=0-S8 
from the energy distribution of the B3V component and (iii) Eb-v=0-87 from the progression of extinction 
along the line of sight. The adopted Eb-v =0.87 ±0.01 is also the amount required by fitting the progenitor with 
theoretical isochrones of appropriate metallicity. The distance is estimated from (a) the galactic kinematics of 
the three components of the interstellar lines, (b) the amount of extinction vs the HI column density and vs 
the dust emission through the whole Galaxy in that direction, from (c) spectrophotometric parallax to the B3V 
companion, from (d) comparison of the observed color-magnitude diagram of field stars with 3D stellar population 
models of the Galaxy, from (e) comparison of theoretical isochrones with the components of the binary system 
in quiescence and found to be around 10 kpc. Pre-outburst optical and IR energy distributions show that the 
component erupting in 2002 was brighter and hotter than the B3V companion. The best fit is obtained for a 
50 000 K source, 0.5 mag brighter than the B3V companion. The latter passed unaffected through the outburst, 
which implies an orbital separation wide enough to avoid mass exchange during the evolution of the binary system, 
and to allow a safe comparison with theoretical isochrones for single stars. Such a comparison suggests that the 
progenitor of the outbursting component had an initial mass ~65 Mq , that it was approaching the Carbon ignition 
stage in its core at the time it erupted in 2002 and that the age of the V838 Mon binary system is close to 4 
million yr. The 2002 event is probably just a shell thermonuclear event in the outer envelope of the star. 

Key words. Stars: evolution - Stars: early type - Stars: individual: V838 Mon - Stars: winds, outflows - ISM: 
dust, extinction - ISM: kinematics 



1. Introduction 

V838 Mon rose from obscurity in early January 2002, 
when it was discovered in outburst by Brown (2002). The 
unusually cool spectrum and strange light-curve helped to 
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keep attention focused on the object for the next three 
months, until the discovery in late March by Henden, 
Munari and Schwartz (2002) of a light-echo rapidly de- 
veloping around V838 Mon. The presence of the first 
Galactic light echo in ~70 years fostered a massive, multi- 
wavelength observing campaign for V838 Mon. A high 
spatial resolution imaging series of the light-echo ex- 
pansion and evolution was collected with HST by Bond 
et al. (2003), recently expanded by new images secured 
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Fig. 1. Optical and IR light-curves of the outburst of 
V838 Mon as recorded with USNO 1.0m and 1.55m tele- 
scopes; some early-time observations come from other 
sources as identified in Munari et al. (2002b). The out- 
burst onset is around Jan 1, 2002, and the first and last 
observations in the picture are for the dates Jan 9, 2002 
and May 11, 2004, respectively. The gaps in the opti- 
cal light-curves correspond to seasonal invisibility. The 
dottcd-dashed lines give the values for the quiescence prior 
to the outburst. 



within the Hubble Heritage Program. An account of the 
spectroscopic, photometric and polarimetric evolution of 
V838 Mon during the first season of visibility was pre- 
sented by Munari et al. (2002a). A major observational 
constraint was the discovery by Desidera & Munari (2002) 
and Munari et al. (2002b) that V838 Mon is a binary sys- 
tem containing a normal B3V star, implying a young age 
and M>7 M Q for its outbursting companion. Over the 
last two years, V838 Mon has been the topic of numer- 
ous refereed papers and conference reports, with a sample 
given below. IR photometry and spectroscopy was pre- 



Table 1. Open clusters within 5° of V838 Mon from the 
Lynga (1987) catalog. 9 is the angular distance in degrees 
from V838 Mon. Reddening and distances are taken from 
the Dias et al. (2002) compilation. 







cluster 


lung. 


gctl. lflt. 


() 


dist. 


lop" ap"e 






4 




(°) 


\ ) 


(°) 


(kpc) 








6 


Bochum 3 


218.79 


+00.35 


1.21 


1.76 


7.890 


0.240 






NGC 2311 


217.76 


-00.69 


1.75 


2.29 


8.600 


0.330 






NGC 2338 


219.88 


+01.01 


2.09 










8 


Haffner 3 


219.82 


-00.01 


2.29 












Berkeley 34 


214.18 


+01.88 


3.71 










10 


NGC 2309 


219.84 


-02.25 


3.88 


2.51 


8.400 


0.350 






NGC 2306 


219.69 


-02.60 


4.12 










12 


NGC 2286 


215.31 


-02.27 


4.16 


2.90 


8.883 


0.170 






Berkeley 77 


219.36 


+05.17 


4.40 








f 




NGC 2302 


219.30 


-03.12 


4.43 


1.18 


7.847 


0.207 






NGC 2323 


221.67 


-01.33 


4.55 


0.93 


8.096 


0.213 




8 


NGC 2324 


213.45 


+03.30 


4.89 


3.81 


8.630 


0.127 






Berkeley 37 


217.23 


+05.93 


4.92 









0.4 



0.3 



a 
p 

<D 
XI 



0.2 



0.1 



12 3 4 

distance (kpc) 

Fig. 2. Extinction as function of distance for open clusters 
within 5° from V838 Mon. 



sented and discussed by Banerjee & Ashok (2002), Crause 
et al. (2003), Evans et al. (2003) and Rushton et al. (2003), 
optical polarimetry and spectropolarimetry by Wisniewski 
et al. (2003a, b) and Desidera et al. (2004), optical pho- 
tometry by Kimeswenger et al. (2002), Barsukova et al. 

(2002) , Kato (2003) and Goranskii et al. (2004), and op- 
tical spectroscopy by Munari et al. (2002a, c), Goranskii 
et al. (2002), Kolev et al. (2002), Osiwala et al. (2002), 
Kipper et al. (2004). Modeling of the outburst was pre- 
sented by Retter & Marom (2003), Soker & Tylenda 

(2003) and Boschi & Munari (2004), while modeling of 
the light-echo expansion was performed by Bond et al. 
(2003), Sugerman (2003) and Tylenda (2004). 

Common to many of these papers were the issues 
of reddening, distance and nature of the progenitor of 
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V838 Mon. However, none of these papers focused specif- 
ically on these topics, and limited their discussion to the 
derivation or adoption of approximate values based on as- 
sumed energy distribution of the outbursting component. 
The aim of the present paper is to attempt an accurate 
derivation of reddening and distance of V838 Mon by care- 
ful comparison of different techniques and approaches that 
do not involve the outbursting component, as well as to 
derive by comparison with theoretical models the nature 
and evolutionary state of V838 Mon prior to the onset of 
the outburst. 

Optical and IR photometry used in this paper was 
obtained with CCD-equipped USNO Flagstaff Station 
(NOFS) telescopes. The UBVR c Ic data were obtained 
with the 1.0m telescope, the JHK'L' with the 1.55m 
telescope. The optical photometry is strictly tied to the 
Landolt (1983, 1992) system of equatorial standards, 
while the infrared photometry uses standard JHK' filters 
and differential measures with respect to local 2MASS 
stars. The low resolution optical spectrophotometry of 
V838 Mon was obtained with the AFOSC+CCD im- 
ager +spectrograph of the 1.82m telescope operated in 
Asiago by the Astronomical Observatory of Padova. At 
the same telescope we secured high resolution spectra 
with the Echelle+CCD spectrograph. Other high resolu- 
tion spectra of V838 Mon were obtained with the Coude 
spectrograph+CCD of the Rozhen 2m telescope, with 
GIRAFFE at the 1.93m SAAO telescope, with SARG at 
the 3.5m TNG telescope, and with FEROS at the 2.2m 
ESO telescope in La Silla. Further details will be provided 
below where necessary. 

Figure 1 presents a review of the optical and IR photo- 
metric evolution of V838 Mon, updated to the end of the 
third season of visibility in May 2004. With the exception 
of some early-time datapoints, the light-curves are entirely 
based on data from the USNO telescopes and therefore 
they arc highly consistent and free from systematic effects 
of color transformation between different local photomet- 
ric systems that would badly affect observations from dif- 
ferent observatories of an object with such extreme colors 
as V838 Mon. 

2. Evidence of patchy extinction 

V838 Mon lies close to galactic equator (b=+l°) and in 
the anti-center quadrant (7=218°). I spite the line of sight 
crosses the Orion, Perseus and Outer spiral arms, inspec- 
tion of the Palomar plates, of the 2MASS, IRAS, UKST- 
Ha maps and the Neckel and Klarc (1980) extinction 
charts do support an apparently smooth star counts dis- 
tribution for several degrees around V838 Mon position. 
Therefore in principle, open clusters and field stars with 
accurate photometry and spectroscopic classification can 
be used to search for a relationship between reddening and 
distance in the direction of V838 Mon. 

Within 5° from V838 Mon there are 13 open clus- 
ters in the catalog of Lynga, (1987). They are listed in 
Table 1. According to the compilation of Dias et al. (2002), 
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Fig. 3. The HD stars closer than 3° from V838 Mon 
are plotted. The crosses mark those within 0.6° from 
V838 Mon. Only HD stars classified by the Michigan 
Project and with Tycho-2 Bt,Vt data are considered. 

only half of them have published determinations of dis- 
tance and reddening. The placing of these 7 clusters on 
a distance-reddening diagram (cf. Figure 2) shows a large 
scatter and basically no obvious trend over such angular 
distances. 

The best source of homogeneous spectroscopic classi- 
fication of field stars is the Michigan Project, which is a 
continuing program at the University of Michigan where 
all HD stars are being reclassified on the MK system, start- 
ing at the south Galactic pole, and ending at the north 
Galactic pole. Volume 5 (Houk & Swift 1999) covers the 
declination strip —12° < S < 5°, thus including the region 
around V838 Mon. 

We selected all HD stars within a radius of 3° from 
V838 Mon and retrieved their Bt and Vr magnitudes 
from the Tycho-2 catalog (H0g et al. 2000). The latter 
were transformed into corresponding Johnson Bj and Vj 
following Besscll (2000) transformations. The reddening of 
each star was than derived by comparison with Fitzgerald 
(1970) intrinsic colors, and the distance computed using 
the absolute magnitudes of the Michigan Project scale 
(from N. Houk web page 1 ). 

Figure 3 plots the distribution of HD field stars in a 
distance-reddening diagram. No well defined progression 
of the reddening with distance is evident in the figure. The 
amount of scatter on the diagram is much larger than jus- 
tifiable in terms of measurement and classification errors, 
and it clearly confirms the evidence from open clusters 
that the reddening is very patchy over the sky area cen- 

1 http:/ /www.astro. lsa.umich.edu/users/hdproj/ 
mosaicinfo/absmag.html 
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Table 2. Absorption along the line of sight to V838 Mon 
according to the analysis of the color magnitude diagram 
of the field stars presented in Figures 4 and 5. 



dist (kpc) 0.1 1.0 2.5 3.0 4.0 5.0 7.0 10.0 
A v 0.01 0.15 0.3 0.7 1.4 1.6 2.1 2.7 



tered on V838 Mon in spite of the smooth star counts 
distribution over the area. Even restricting to the eleven 
HD stars closer than 0.6° (marked by crosses in Figure 3) 
does not help to define a well-behaved relation between 
distance and reddening. This indicates that the redden- 
ing coherence area around V838 Mon has a significantly 
smaller radius. 

3. The field stars close to V838 Mon 

Noting the patchy reddening structure over degree-sized 
areas, we have instead focused on a much smaller angular 
scale in looking for a coherent relation between distance 
and reddening along the line of sight to V838 Mon. Deep 
UBVRqIg photometry of field stars within a 11x11 ar- 
cmin field centered on V838 Mon has been obtained with 
the NOFS 1.0m telescope. The observed color-magnitude 
diagrams of these field stars have then been compared to 
and fitted with synthetic diagrams from stellar population 
models of our Galaxy. The photometric data of the field 
stars are electronically available 2 . 

V838 Mon lies close to the galactic equator and in the 
general anti-center galactic direction, where stellar counts 
are largely dominated by the disk population. Bertelli et 
al. (1995) showed that the slope of the main sequence of 
the disk population in color-magnitude diagrams is mainly 
governed by the extinction along the line of sight. At any 
V magnitude the bluest stars on the main sequence can be 
interpreted as the envelope of the main sequence turnoffs 
of the population having absolute magnitude M tur , shifted 
toward fainter magnitudes and redder colors by the in- 
creasing distance and corresponding extinction. Starting 
from an initial guess, the amount of extinction at increas- 
ing distances is adjusted (via x 2 test on the color distribu- 
tion as function of magnitude) until a satisfactory agree- 
ment between the main sequence blue edge location in the 
observations and in the model is reached. The distribution 
of observational errors is evaluated and taken into account 
in the process. 

The Galaxy model that we used is based on the code 
described by Bertelli et al. (1995) and recently revised by 
Vallenari et al. (2000, 2004). The stars are generated ac- 
cording to the evolutionary model of the Galaxy, and then 
are distributed along the line of sight following the spa- 
tial model of the Galaxy. All stellar populations are taken 
into account. The generation of the synthetic population 
makes use of the set of the stellar tracks by Girardi et 

2 At http://ulisse.pd.astro.it/V838_Mon/ 
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Fig. 5. Extinction along the line of sight to V838 Mon 
(thick and solid line) as implied by the galactic synthetic 
population fit in Figure 4 to the observed color-magnitude 
diagram of field stars. The dotted line represents the 
Mendez & van Altena (1998) model for the dust in the 
Galaxy, the short dashed-dotted line is the Drimrnel et al. 
(2003) model, while the dashed lines give the total galactic 
extinction along the line of sight according to the Schlegel 
et al (1998) COBE maps. The long dashed-dotted line is 
the extinction behavior according to the Neckel & Klare 
(1980) maps (poorly populated at >3 kpc). 

al. (2000). A double exponential mass distribution inside 
the thin disk is adopted having a scale length of 2500 pc 
and a scale height varying with the age of the population, 
ending at 350 pc for the oldest component. A thick disk 
component is included with a scale height of 800 pc. A 
stochastic age metallicity relation for the disk is assumed, 
with Z going from 0.008 to 0.03. 

We performed the analysis in both V,B-V and V, V-I 
planes, with similar results. Figure 4 presents the distri- 
bution of observed and synthetic stars in the V,V-I dia- 
gram, and Figure 5 gives the resulting distribution of the 
extinction (Ay) along the line of sight. We experimented 
with different star formation histories in both the thin 
and thick galactic disks, but Ay=f(d) never changed by 
more than 0.1 mag at any point along the line of sight. 
The values reported in Table 4 and adopted here are those 
providing the best results in the % 2 tests. The error bars 
in Figure 5 include both the effect of observational er- 
rors and incompleteness as well as different star forma- 
tion histories in the galactic disk over the last 10 Gyr. 
The magnitude of turnoff stars beyond 6 kpc is at V >19, 
where observational errors become significant and com- 
pleteness begins to be affected. Therefore, the results in 




Fig. 4. The left panel plots on a V, V — I diagram the stars within a 11x11 arcmin field centered on V838 Mon 
as observed with the NOFS 1.0m telescope. The right panel represents the fit to the observed distribution using the 
Padova Galaxy Model, normalized to the same number of stars (1358). 



Table 4 and Figure 5 are less accurate beyond 6 kpc. To 
increase the accuracy at larger distances it is advisable 
that much deeper photometry is obtained in future ob- 
servations with 8m class telescopes. Comparing the HR 
diagrams of field stars in different colors it is clear that a 
standard Ry = Ay/Es-v =3.1 extinction law applies. 



In Figure 5 comparison is made with other Ay=f(d) 
already presented in literature. Mendez & van Altena 
(1998) make use of the large-scale properties of the dust 
layer in the Galaxy to derive the absorption in the plane of 
the Galaxy. Their model indicates that a differential opti- 
cal absorption of 0.5 mag kpc -1 is adequate to reproduce 
the available reddening maps in the range 2<r<6 kpc. The 
internal uncertainties reported by Mendez & van Altena 
are of the order of i?B_y=0.25, or Ay— 0-77 mag for a 
standard Ay —3.1 xEb-v- Schlegel et al (1998), using the 
COBE-DIRBE 100 and 240 fim all sky survey, derived a 
map of dust emissivity, which was converted into a map 
proportional to dust column density and, finally, to a map 
of the total amount of reddening through any direction in 
the Galaxy. These maps are expected to suffer in accu- 
racy close to the galactic plane, at |6| < 5°. Drimmel et 
al (2003) recently presented a dimensional model of the 
dust distribution based on COBE-DIRBE infrared data 
which aims to be useful closer to the Galactic plane, since 
it accounts for the presence of spiral arms. As stated by 
Drimmel et al., regions having anomalous emission due to 
warm dust (like Orion or the Galactic center) are not well 
described by their model. Neckel & Klare (1980) extinc- 
tion maps are based on observed B — V colors of stars 
of known spectral type and luminosity class so that a 
spectro-photometric parallax can be estimated. Neckel & 
Klare data are only tentative at distances larger than 3 kpc 
given the paucity of data. 



4. The Nal and Kl interstellar lines 

The interstellar lines are a powerful mean for deriving the 
reddening along the line of sight. Munari & Zwitter (1997) 
have calibrated accurate relations between reddening and 
equivalent width of Nal and KI interstellar lines that we 
apply to V838 Mon in this section. 

We have secured a number of high resolution spec- 
tra during the outburst of V838 Mon which are suitable 
to reveal and measure the interstellar Nal and KI lines. 
They arc summarized in Table 3 (available in electronic 
form only). The general appearance of interstellar lines 
in V838 Mon spectra is illustrated in Figure 6. There are 
three main components contributing to the profiles of in- 
terstellar lines. They are seen well separated in the KI 
line observed at a resolving power 48 000 with FEROS at 
the 2.2m at ESO. The two components at +29 and +42 
km sec -1 are blended at the lower resolution of the other 
instruments, particularly in the case of the stronger Nal 
lines. For this reason these lower-resolution spectra show a 
narrow red component, at +65.1 km sec" 1 radial velocity, 
and a broader blue component, at +37.7 km sec" 1 radial 
velocity. It is worth noting that the average velocity of the 
two blue KI components (at +29 and +42 km sec -1 ) in 
the FEROS spectrum, weighted according to their equiv- 
alent widths, matches exactly the velocity of the blended 
blue component (+37.7 km sec -1 ) in the lower resolution 
spectra of both Nal and KI lines. 

The results of the measurement of heliocentric wave- 
length, equivalent width and FWHM of the Nal and KI 
lines of spectra listed in Table 3 are given in Table 4 (avail- 
able in electronic form only), separately for the blue and 
the red components. Their mean values and standard de- 
viation are given in Table 5. Extensive tests on the spec- 
tra suggests that the minimal apparent changes with time 
traceable in Table 4 are not the result of intrinsic vari- 
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Fig. 6. Interstellar lines in the high resolution optical 
spectra of V838 Mon. The upper panel shows the profile of 
the Nal Dl line superimposed on the broad P-Cyg stellar 
line. The profile is composed of a narrow red component 
and a blue unresolved blend. The blue blend is resolved in 
a higher resolution FEROS observations of the KI 7699 A 
line of V838 Mon on the lower panel. 

Table 3. Journal of the high-resolution spectroscopic 
observations with the Asiago, SAAO, ESO, TNG and 
Rozhen telescopes (available in electronic form only). 



Table 4. Results of the multi-Gaussian fitting to the Nal 
and KI interstellar lines in the high resolution spectra of 
V838 Mon listed in Table 1 (available in electronic form 
only). 



ability of the interstellar lines themselves. Instead, they 
arise from the difficulty (arbitrary) in tracing the level of 
the underlying continuum in the continuously changing 
shape of the underlying broad stellar P-Cyg profiles (cf. 
Figure 6), as well as the large assortment of spectrographs 
- each with its own PSF - contributing to the monitoring 
of the interstellar lines. 

The time behavior of heliocentric radial velocity equiv- 
alent width and FWHM is plotted in Figure 7. It clearly 
indicates that since earliest outburst phases these val- 
ues have remained constant, being unaffected by the 
large changes in the mass loss (amount and velocity) of 
V838 Mon during its active state. This indicates that the 
ejected material had not reached and swept away any pre- 
existing circumstellar material giving rise to Nal and KI 
absorptions, which is consistent with the absence of X- 
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Fig. 7. Heliocentric radial velocities, equivalent widths 
and full width at half maximum from Table 4 of the inter- 
stellar lines as function of time since the onset of the out- 
burst of V838 Mon around Jan 1, 2002. Their constancy 
suggests that the violent eruption has not interfered with 
possible pre-existing circumstellar material. 

ray detection in Chandra observation of V838 Mon by 
Orio et al. 2003 and the central dust-free region described 
by Tylenda (2004). All three components of the interstel- 
lar lines seem truly interstellar and not circumstellar be- 
cause (a) ground-based and HST images of the light-echo 
(Munari et al. 2002c, Bond et al. 2003, and a recent Hubble 
Heritage Program image) show a hole surrounding the 
central star, (b) the reddening from interstellar lines (cf 
next section) matches that from field stars, and above all 
because (c) the radial velocities of the three components 
of the interstellar lines (65, 42 and 29 km sec -1 ) precisely 
correspond to those of the HI complexes (68, 43 and 27 
km sec -1 ) observed along the line of sight to V838 Mon 
in the Leiden/Dwingcloo HI Survey (Hartmann & Burton 
1997). 

4.1. Reddening from interstellar lines 

From the wavelength of maximum polarization, Desidera 
et al. (2004) have demonstrated how the reddening affect- 
ing V838 Mon follows the normal Ry = Ay/Es-v =3.1 
extinction law. The same result follows from the anal- 
ysis of UBVRI photometry of field stars in sect. 3. We 
can therefore safely apply to interstellar lines the Munari 
& Zwitter (1997) calibration for the average interstel- 
lar medium to the equivalent width of Nal and KI lines 
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Table 5. Mean values and standard deviations of the ra- 
dial velocity (RV), full width at half maximum (T) and 
equivalent width (E.W.) of the Nal D and KI interstellar 
lines in the spectra of V838 Mon detailed in Table 2. 
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Nal Dl 


+38.4 
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29.8 1.7 


0.610 


0.035 


Nal D2 


+36.4 


1.8 


29.0 1.1 


0.625 


0.044 


KI 


+38.3 


2.1 


18.1 3.4 


0.100 


0.054 








red component 






Nal Dl 


+65.7 


1.8 


18.3 1.7 


0.329 


0.038 


Nal D2 


+64.3 


1.5 


19.4 0.9 


0.366 


0.009 


KI 


+65.4 


1.8 


13.4 3.4 


0.128 


0.017 



in Table 5. The resulting reddening is Eb-v=0-86 from 
KI and Eb-v—0-&5 from Nal. The results are in excel- 
lent agreement between the two sets of lines. We adopt 
Eb-v=0-86 as the reddening indicated by the interstel- 
lar lines, giving higher weight to the results from KI lines 
which are still far from saturation compared to the nearly 
saturated Nal lines. 

4.2. Distance from radial velocity of interstellar lines 
and the galactic rotation 

Some constraints on the distance to V838 Mon can be 
derived by combining the observed radial velocity of in- 
terstellar material and the rotation curve of the Galaxy. 

The three principal components of the interstellar ab- 
sorption lines have heliocentric radial velocities of +65.1, 
+41.9 and +28.7 km sec -1 . In Figure 8 they are over- 
plotted onto the rotation curve of the Galaxy from Hron 
(1987) and onto the velocity field in the direction of 
V838 Mon as derived by Brand & Blitz (1993). 

The velocity of the two components at +41.9 and 
+28.7 km sec -1 is compatible with a broad range of dis- 
tances, while the component at +65.1 km sec" 1 supports a 
distance to the absorbing medium not shorter than 3 kpc, 
with a corresponding larger distance for V838 Mon itself. 
Similar lower limits to the distance to V838 Mon were de- 
rived by Wisniewski et al. (2003a) and Kipper et al. (2004) 
from similar kinematic arguments based on the observed 
velocity of interstellar lines. 

5. The B3V companion 

Munari et al. (2002c) commented on the fact that the 
UBVRcIc light-curves of V838 Mon started to exhibit a 
puzzling bluing in U — B 90 days and in B — V 115 days 
after the outburst onset, while V — I was rising to very 
red values, approaching those of brown dwarfs. The reason 
became apparent when Desidera and Munari (2002) dis- 
covered on spectra at day 274 that the contribution from 
the very cool outbursting component was less in the blue 
part of the spectrum, revealing the presence of a fainter 



and hotter companion. V838 Mon was thus shown to be 
a binary system, and the hot companion was classified 
by Munari et al. (2002b) on spectra for day 300 as a B3 V 
star. The presence of the hot companion was confirmed on 
higher resolution spectra by Wagner and Starrficld (2002). 
The composite nature of optical spectra of V838 Mon dur- 
ing this phase is illustrated in the top panel of Figure 9. 
Later on there was an increase in the contribution of the 
outbursting component at shorter wavelengths, as illus- 
trated in the bottom panel of Figure 9. 

The average color of V838 Mon at the time of the 
naked visibility of the B3V companion was B V=+0.68 
(cf. Figure 1 and Munari et al. 2002b). Comparing with 
the intrinsic colors of a B3V from Fitzgerald (1970) it 
results in £^B_y=0.88, in excellent agreement with the 
results from interstellar Nal and KI lines. 

The metallicity in the galactic disk at the galactocen- 
tric distance of V838 Mon is [Fe/H]=-0.7, about half dex 
lower than in the solar neighborhood (e.g. Friel et al. 2002, 
Hou et al. 2003). The effect on the B-V color of the B3 V 
component in V838 Mon of such a reduction in metallic- 
ity is a very minor one, being on the Raylcigh- Jeans tail 
of the energy distribution of a hot star. Integrating the 
transmission profiles of the Landolt's B and V bands (to 
which our photometry is tied) to the 2500-10500 A syn- 
thetic Kurucz spectral library of Munari et al. (2004), the 
net effect is just 0.007 mag. 

The average V magnitude of V838 Mon at the 
time the V band was dominated by the radiation 
from the B3V component is ^=16.05+0.05. Coupled 
with the _Eb_v=0.87±0.01 and an absolute magnitude 
M y =-1. 70+0.05 for a B3V star from Houk (2004), it 
implies a distance of 10 kpc to V838 Mon. 



6. Combined results 

6.1. Reddening 

Two assumption-independent methods, the interstellar 
atomic absorption lines and the colors of the B3V com- 
panion, provide consistent results on the reddening affect- 
ing V838 Mon: £ B _y=0.86 and £ B _ V =0.88, respectively. 
There is no differential extinction between the outbursting 
star and the B3V companion. 

The same amount of reddening is derived from the 
modeling of the HR diagram of field stars close to 
V838 Mon performed in sect. 3. In fact, entering Table 2 
or Figure 5 with the distance (10 kpc) to V838 Mon from 
the spectro-photomctric parallax to the B3V companion, 
i?B-v=0.87 is obtained. 

We therefore conclude that the extinction toward 
V838 Mon follows the standard R Y = A v /E B -v =3.1 
law and the reddening amounts to E B -v =0.87+0. 01. 
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Fig. 8. The radial velocities of the three components 
of Nal and KI interstellar lines (+28.7, +41.9 and 65.1 
km sec -1 ) are over-plotted on the velocity field from 
Brand & Blitz (1993) and on the rotation curve of the 
Galaxy of Hron (1987). For the latter two additional 
curves are plotted, corresponding to the extremes of the 
allowed ranges for A and a parameters. 

6.2. Distance 

There are several and independent arguments and evi- 
dences that support a long distance scale to V838 Mon, of 
the order of 10 kpc. 

Analysis of the high-resolution HST polarimctry im- 
ages of the light-echo led Bond et al. (2003) to place a 
lower limit of 6 kpc to the distance of V838 Mon. Working 
on the same HST material, Tylenda (2004) revised the dis- 
tance to 8±2 kpc. 

The spectro-photometric parallax to the B3V compan- 
ion derived in sect. 5 is 10 kpc. Furthermore, reversing the 
argument of the previous section, we can enter Table 2 or 
Figure 5 with the reddening £?b-v=0.87 (corresponding 
to Av=2.7) determined from the interstellar lines and the 
colors of the B3V companion, obtaining the same distance 
of 10 kpc to V838 Mon. 

Two other arguments are in favor of a large dis- 
tance to V838 Mon, given its proximity to the galactic 
plane (b=— 3°51' galactic latitude). First, we have seen 
above that the HI radio observations in the direction of 
V838 Mon reveal three components whose velocity match 
exactly that of the three components seen in the interstel- 
lar absorption lines. Given the fact that the radio obser- 
vations integrate along the whole line of sight through the 
Galaxy and that no HI is seen beyond V838 Mon (in the 
sense that no corresponding interstellar lines are detected 
in the high- resolution spectra), it is straightforward to 
conclude that V838 Mon lies at great distance with most of 
the galactic gas, in that direction, in front of it. Second, 
Schlcgel et al. (1998) extinction maps are based on the 
amount of infrared emission by the dust integrated along 
the whole line of sight through the Galaxy. The amount of 
extinction affecting V838 Mon (Eb-v=0-87 correspond- 
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4000 5000 6000 7000 

wavelength (A) 

Fig. 9. Comparison between Asiago AFOSC spectra of 
V838 Mon taken in October 2002 and December 2003. At 
the earlier date the outbursting component is only visible 
toward the red, revealing the presence of a B3 V compan- 
ion whose emission dominates the U, B and V bands at 
that time. The later spectrum illustrates how the emission 
of the outbursting component has retraced back to shorter 
wavelengths during the 2003 /04 visibility season, a change 
already described in its early development by Wagner et 
al. (2003). 

ing to Ay=2.7 mag) match the value of Schegel et al. 
maps in that direction, supporting again the notion that 
V838 Mon lies at large galactocentric distances. 

All these independent determinations support the con- 
clusion that V838 Mon lies in the outer part of the disk of 
the Galaxy, at a distance of ~10 kpc from the Sun, corre- 
sponding to a galacto-centric distance of <~17.5 kpc and a 
height above the galactic plane of ^650 pc. 

Tylenda (2004) labeled as "naive interpretation" the 
early distance estimates of Munari et al. (2002a) and 
Kimeswcnger et al. (2002) based on the first determina- 
tions of the angular expansion rate of the light-echo on 
early ground-based discovery images. The Munari et al. 
(2002a) implicit (but quite obvious) assumption was that 
the light-echo was originating in a circumstellar disk seen 
pole-on. Such an assumption was based on the fact that 
the Nal and KI lines were not tracing a circumstellar com- 
ponent, something to be expected in the case of a ho- 
mogeneous spherical distribution of material centered on 
the object itself. It was only much later that high res- 
olution HST imaging revealed that there is a clear void 
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Fig. 10. The crosses represent the energy distribution 
of V838 Mon at the time of naked visibility of the hot 
companion (V=16.05, B-V=+0.68, £/-£=-0.06) com- 
pared with a Kurucz's spectrum for a metal poor B3V 
star (T cff =19000 K, lg 5=4.0, [Fe/H]=-0.7) reddened by 
£?B_y=+0.87. The other symbols give the energy distri- 
bution of V838 Mon in quiescence: stars from Goranskii 
et al. (2004), triangles from Kimeswenger et al. (2002), 
squares from DENIS survey, dots from 2MASS survey, 
circles from USNO-B magnitudes re-calibrated against 
the comparison sequence of Munari et al. (2002), giving 
i?=15.28 and i?c=14.22. The quiescence energy distribu- 
tion is fitted with the combination of Kurucz's spectra for 
the B3V stars and the hotter companion. 



of circumstcllar material precisely along the line of sight 
to V838 Mon, regardless of the true 3D shape of the cir- 
cumstellar dust giving rise to the light-echo. So, the now 
"naive" approach was reasonable at the time. 

6.3. Progenitor 

Figure 10 shows the energy distribution of the B3 V com- 
panion in the UBV bands (<V>=16.052, <5>=16.736, 
<£/>=16.676 average values from USNO 1.0m photom- 
etry described in sect.l) compared with a Kurucz's syn- 
thetic spectrum (from the library of Munari et al. 2004) 
with parameters appropriate for the B3V companion 
(T eff =19 000 K, log .9=4.0, [Fe/H]=-0.7). The match is ex- 
cellent. In the same figure, the pre-outburst brightness of 
V838 Mon from various sources are plotted. In the infrared 
they come from the 2MASS and DENIS surveys, while the 
B and Rc values are estimates of the same POSS-I and 
SERC plates according to different authors that used dif- 
ferent calibrations. The entries at / band (note the slightly 
different effective wavelengths) come from the DENIS and 
POSS-II surveys. 

The figure clearly shows that the progenitor of the out- 
bursting component was brighter and hotter than the B3 V 
companion. By exactly how much it depends on the still 
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Fig. 11. Padova isochrones (from Bressan et al. 1993, 
Fagotto et al. 1994, Girardi et al. 2000) for Z=0.004 and 
three ages scaled to the distance (10 kpc) and reddening 
(£ B -v=0.87) of V838 Mon. The isochrones include the 
effect of mass loss and are corrected for the reddening- 
induced distortion described in Fiorucci & Munari (2003, 
their Figure 8). The open circle indicates the position of 
the V838 Mon component that erupted in 2002 if its tem- 
perature in quiescence was 50 000 K, and the arrow gives 
the shift if the temperature was 25 000 K (cf. the energy 
distribution fit in Figure 10). 

uncertain magnitude of V838 Mon on pre-outburst photo- 
graphic plates. While the brightness seems well established 
in the infrared by the excellent agreement of DENIS and 
2MASS data, it presents a large scatter in the estimates 
by different authors of the same POSS and SERC plates. 
It is clear from Figure 10 that a careful re-evaluation of 
all existing pre-outburst photographic material would be 
beneficial in establishing the nature of V838 Mon prior 
to the outburst. Such a detailed analysis must first pro- 
vide a perfect comparison sequence, a tight control over 
the plate measurement and its modeling, a careful knowl- 
edge of the exact band-pass profile of the given emulsion- 
filtcr-telescope-atmosphere combination and plate devel- 
oping conditions, as well as careful application of color 
equation corrections. 

The best fit to the V838 Mon quiescence data in 
Figure 10 is the B3V plus a 50000 K star with V==15.55 
and B — V=+0.535. Both stars are equally reddened by 
Eb-v=0-87. The combined spectrum of the two stars 
is plotted in Figure 10 as a solid line. Ignoring our re- 
estimate of the magnitude of V838 Mon on the POSS and 
SERC plates listed in the caption of Figure 10, a some- 
what lower temperature for the eruption progenitor can 
be obtained, with a lower limit of 25 000 K. The combina- 
tion of the spectra of the B3V plus a 25 000 K star with 
V = 15.74 and B - V=+0.634 is plotted as a dotted line in 
Figure 10. 

The location of the components of the binary sys- 
tem V838 Mon is compared with theoretical isochrones 
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in Figure 11, for the [Pe/H]=— 0.7 metallicity appropriate 
to the galactocentric distance of V838 Mon and for the so- 
lar case as a matter of reference. The Padova theoretical 
isochrones have been scaled to a distance of 10 kpc and 
E B _ v —0.87 reddening following the standard i?y=3.1 ex- 
tinction law. The reddening dependent deformation over 
the HR plane discussed and quantified by Fiorucci & 
Munari (2003, cf. their Figure 8) for the UBVR C I C pho- 
tometric system (Vilnius and USA reconstructions) has 
been applied to the isochrones in Figure 11. 

As displayed by Figure 11, there is a perfect match of 
the [Fe/H]=— 0.7 isochrones with the observed magnitude 
and colors of the B3 V component, with the progenitor of 
the outbursting companion lying on the isochrone for an 
age of 4 million yr, close to the position of the central 
ignition of Carbon for an initial mass of ~65 M Q . The 
progenitor lies close to the region where Wolf-Rayet stars 
are usually located. It would be harder to fit with theo- 
retical evolutionary tracks the position of the progenitor 
of the outbursting component if its temperature would be 
decreased to the lower limit of 25 000 K of Figure 10, in- 
dicated by the head of the arrow in Figure 11. 

We therefore conclude that the progenitor of the out- 
bursting component of the binary system V838 Mon was 
hotter, brighter and born far more massive than the 7 M 
B3V companion, of an age not far from 4 million yr. In 
this respect, the circumstellar scattering material giving 
rise to the light-echo is probably the result of the mas- 
sive mass loss that such massive objects experience, where 
^half of the mass is lost during the main sequence phase. 
Given the location of V838 Mon in the outskirts of the 
Galactic disk at galactic longitude Z=218°, the large mass 
and young age could sound problematic. However, young 
clusters and massive stars at great distances in the anti- 
center direction are already known. For ex., Fitzgerald & 
Moffat (1976) reported a distance of ^7 kpc for the very 
young cluster Ruprecht 44 (Z=245°) rich in O and WR 
stars, or Marco et al. (2001) derived a distance of 6 kpc 
and an age of 4 Myr for the cluster NGC 1893 (1=173°) 
that harbors 05V member stars (which mass is ~65 M Q 
according to Straizys & Kuriliene (1981) tabulation). 

The outburst experienced in 2002 does not appear as 
the terminal event in the life of the massive progenitor, 
but instead more probably as a thermonuclear shell flash 
in the outer layers of the star as could be expected in the 
case of He after most of the H-rich outer envelope has been 
blown away by the strong wind that characterize this type 
of stars. A detailed analysis of the evolutionary status of 
the progenitor will be presented elsewhere (Munari et al., 
in preparation). 
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Table 3. Journal of high-resolution spectroscopic obser- 
vations. HUT = heliocentric UT date (dd/mm/yy); HJD 
— corresponding heliocentric JD (—2452000); N= num- 
ber of spectra; Exp^= total exposure time (in sec); S/N 
= signal-to-noise of the continuum around Nal and/or KI 
in the combined spectrum; D — dispersion; Res = resolv- 
ing power (A/AA); slit — aperture on the sky (in arcsec) 
of the spectrograph entrance slit; lines = interstellar lines 
observed. 
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363.28 


1 


1800 


205, 230 


0.171, 0.225 


20 000 


2.0 


Nal, KI 


Asiago 


22.82/04/02 


387.32 


1 


1200 


23, 70 


0.171, 0.225 


20 000 


2.0 


Nal, KI 


Asiago 


25.71/04/02 


390.21 


1 


1200 


10 


0.075 


39 000 


2.0 


Nal 


SAAO 


22.27/10/02 


569.77 


1 


1600 


14, 23 


0.043, 0.055 


29 000 


1.6 


Nal, KI 


TNG 


19.09/12/02 


627.59 


2 


3600 


15 


0.225 


20 000 


2.0 


KI 


Asiago 


18.94/02/03 


689.44 


2 


3600 


17 


0.225 


20 000 


2.0 


KI 


Asiago 


20.95/02/03 


691.45 


3 


5400 


19 


0.225 


20 000 


2.0 


KI 


Asiago 


25.04/02/03 


695.54 


5 


8100 


11, 33 


0.030, 0.030 


48 000 


1.8 


KI 


ESO 
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Table 4. Results of the multi-gaussian fitting to the Nal and KI interstellar lines in the high resolution spectra 
of V838 Mon listed in Table 1 (rest wavelengths at 5895.923, 5889.953 and 7698.979 A for Nal Dl, D2 and KI, 
respectively). All wavelengths (A) are heliocentric corrected, and the FWHMs (T) are corrected for the width of the 
instrumental PSF. Both are expressed in km sec -1 . The first two digits of the wavelengths (59 for both Nal Dl and D2, 
76 or 77 for KI) are dropped for table compactness. The equivalent widths are in A, and computed relative to a linear 
fit of the stellar continuum adjacent to the lines. The nomenclature blue and red components refers to the broader line 
at <+37.7> and the narrower one at <65.1> km sec -1 , respectively, that compose the profile of all interstellar lines, 
and that are well illustrated by the Nal Dl line at the top of Figure 1 (superimposed on the wide Nal P-Cyg stellar 
profile) . 



Nal Dl blue comp. Nal Dl red comp. Nal D2 blue comp. Nal D2 red comp. KI blue comp. KI red comp. 







A 


E.W. 


r 


A 


E.W. 




r 




A 


E.W. 




r 




A 


E.W. 




r 




A 


E.W, 




r 




A 


E.W. 




r 


25.98/01/02 


96 


.711 


0.642 


0.623 


97.263 


0.365 





.416 


90 


.688 


0.589 





.562 


91 


.241 





,350 





.360 


00 


.057 










00 


.607 










28.91/01/02 


96 


.690 


0.667 


0.664 


97.238 


0.337 





,412 


90 


.694 


0.653 





,600 


91 


.240 





.378 





.402 


99 


.949 










00 


.617 










29.92/01/02 


96 


.702 


0.644 


0.621 


97.234 


0.350 





.399 


90 


.700 


0.703 





,652 


91 


.251 





,355 





.411 


99 


.996 










00 


.625 










25.83/02/02 


96 


.731 


0.608 


0.590 


97.236 


0.289 





,321 


90 


.696 


0.615 





549 


91 


.220 





.376 





.389 


99 


,951 





.086 





.535 


00 


.667 


0, 


,119 


0, 


.343 


27.77/02/02 


96 


.634 


0.602 


0.589 


97.147 


0.289 





,322 


90 


,614 


0.634 





,579 


91 


.151 





.362 





.391 


























1.76/03/02 


96 


.634 


0.611 


0.587 


97.160 


0.304 





.332 


90 


,627 


0.661 





,592 


91 


.167 





.364 





,389 


























1.84/03/02 


96 


.708 


0.650 


0.620 


97.219 


0.285 





,318 












91 


.221 





.362 





.381 


























3.82/03/02 


96 


,703 


0.618 


0.583 


97.220 


0.304 





.322 


90 


,695 


0.676 





,598 


91 


.228 





.367 





.387 


























5.84/03/02 


96 


,698 


0.617 


0.582 


97.218 


0.306 





.321 












91 


.217 





.369 





.375 


























7.86/03/02 


96 


.692 


0.626 


0.589 


97.214 


0.306 





.325 
















































26.79/03/02 


96 


.671 


0.556 


0.535 


97.214 


0.389 





.416 


90 


,654 


0.584 





,534 


91 


.203 





.356 





.372 


99 


,997 





.075 





.567 


00 


.679 


0, 


.107 


0, 


,299 


27.84/03/02 


96 


.682 


0.555 


0.539 


97.229 


0.390 





,425 


90 


.680 


0.605 





.554 


91 


.230 





.379 





.345 


00 


,020 





.064 





.481 


00 


.702 


0, 


.107 


0, 


.320 


29.78/03/02 


96 


.702 


0.618 


0.579 


97.244 


0.328 





,342 


90 


.693 


0.582 





,528 


91 


.242 





.366 





.373 


99 


.847 





.064 





.432 


00 


.567 


0, 


,119 


0, 


,297 


22.82/04/02 


96 


.660 


0.540 


0.539 


97.217 


0.397 





.458 


90 


.606 


0.568 





.510 


91 


.189 





.378 





.383 


99 


,991 










00 


.694 










25.71/04/02 


96 


.640 


0.590 


0.557 


97.207 


0.357 





.331 
















































22.27/10/02 


96 


.596 


0.610 


0.588 


97.168 


0.299 





.327 
























99 


,990 





.149 





.388 


00 


.661 


0, 


,139 


0, 


,300 


19.09/12/02 








































00 


,009 





.202 





.518 


00 


.736 


0, 


,137 


0, 


,364 


18.94/02/03 




















































00 


.622 


0, 


,126 


0, 


411 


20.95/02/03 








































99 


,968 





,061 





.323 


00 


.704 


0, 


,155 


0, 


.530 


25.04/02/03 








































99 


.716 





.039 





.151 


00 


.655 


0, 


,145 


0, 


.240 



00.054 0.065 0.232 



